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I Abstract 

o 

In this paper, we quantify the systematic impact of the non-spherical shape of transiting planets, due to tidal forces and rotation, on the 
! * observed transit depth. Such a departure from sphericity leads to a bias in the derivation of the transit radius from the light curve and 

affects the comparison with planet structure and evolution models which assume spherical symmetry. As the tidally deformed planet 
projects its smallest cross section area during the transit, the measured effective radius is smaller than the one of the unperturbed 
spherical planet (which is the radius predicted by ID evolution models). This effect can be corrected by calculating the theoretical 
shape of the observed planet. 

Using a variational method and a simple polytropic assumption for the gaseous planet structure, we derive simple analytical expres- 
sions for the ellipsoidal shape of a fluid object (star or planet) accounting for both tidal and rotational deformations. We determine the 
characteristic polytropic indexes describing the structures of irradiated close-in planets within the mass range 0.3 Mj < M p < 75 Mj, 
at different ages, by comparing polytropic models with the inner density profiles calculated with the full evolution code. Our cal- 
culations yield a 20% effect on the transit depth, i.e. a 10% decrease of the measured radius, for the extreme case of a lMj planet 
orbiting a Sun-like star at 0.01 AU, and the effect can be larger for smaller mass objects. For the closest planets detected so far (< 0.05 
AU), the effect on the radius is of the order of 1 to 10% (three times more for the mean density), by no means a negligible effect, 
enhancing the puzzling problem of the anomalously large bloated planets. These corrections must thus be taken into account for a 
correct determination of the radius from the transit light curve and when comparing theoretical models with observations. 
Our analytical expressions can be easily used to calculate these corrections, due to the non-spherical shape of the planet, on the 
observed transit depth and thus to derive the planet's real equilibrium radius, the one to be used when comparing models with 
observations. They can also be used to model ellipsoidal variations of the stellar flux now detected in the CoRoT and Kepler light 
1 1 curves. We also derive directly usable analytical expressions for the moment of inertia and the Love number (fo) of a fluid planet as a 

function of its mass and, in case of significant rotation, for its oblateness. 
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1. Introduction These observations motivate us to investigate the deforma- 
tion of the planet with respect to a spherical body, because of 
<~H As me measurement of the radii of close-in transiting planets tidal or rotational forces . while prev ious studies focused on the 
O continues to gain in accuracy, providing stringent constraints det ectability of the oblateness of a flattened body, we address in 
on exoplanet theoretical models, any source of errors in the ra- the present paper the more general problerrii namely the deter . 
r-H dms determination must be determined with precision. Current mination of the genera i shape f „ planet (or star) distorted by 
> ground and space-based photometric observations of the host both a tidal and a cen trifugal potential, and its impact on the tran- 
. ■ rH stars of transiting planets enable us to address new problems. sit depth and thus on the det e rmi nation of its correct radius. In 
^ The first direct detection with Spitzer of the light emitted by order to compute the ellipsoidal shape (flattening and triaxiality) 
H the planet (Deming et al. 2007) opened a new path to probe the of a gaseoufj bodyi we deriye in Sect 2 a simple analytical model 
^ physical properties of the surface and the atmosphere of transit- of the intemal structure of the object (Lai et aL 1994)j based on 
ing exoplanets. Among the first results of the Kepler mission, the polytropic equation of state. The polytropic indexes are cal- 
the detection of ellipsoidal variations of the host star induced by ibrated in Sect 4 by comparing with nume rical models describ- 
tidal interaction with a low mass companion has been claimed ing the stmcture of strong i y ilTa diated gaseous planets (Leconte 
Welsh et al. (2010). More recently, Carter & Winn (2010a,b) et al 2QQ9) In Sect . 3, we present directly usable analytical ex- 
showed that light curve analysis can put direct constraints on the pressions giving the shape (oblateness and triaxiality) of a dis- 
actual shape of transiting planets. They also investigated the im- torted planet (or star) as a function of its mass and po i ytrop i c 
pact of the precession of an oblate object with a non zero obhq- index and compare our estimates with the numerical method 
uity around the orbital axis on the shape and timing of the transit outlined in Appendix A and with the measured values for the 
s ig nai - major planets of our solar system. As a by-product, our model 

yields analytical expression for the first gravitational moment J2 
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and the love number ki of a self gravitating fluid body. Finally, 
Sect. 5 quantifies the effect of the non-sphericity of the planet on 
the transit depth. 

We find that as the planet transits across the stellar disc, we 
only see the smaller cross section of its actual ellipsoidal shape 
so that the depth of the transit is decreased with respect to the 
expected signal for a spherical object, as discussed by Li et al. 
(2010) in the case of WASP-12b. This implies that the radius 
inferred from the light curve analysis, derived under the assump- 
tion of spherical planet and star, underestimates the real equilib- 
rium radius of the object. This bias needs to be corrected for 
a proper comparison with theoretical ID numerical simulations 
of the structure and evolution of extrasolar planets and enhances 
the actual discrepancy between theory and observation for the so 
called "bloated" planets. 

2. Variational method for compressible ellipsoids 

In this section, we briefly describe the energy variational method 
developed by Lai et al. (1993) and Lai et al. (1994) (hereafter 
LRS 1 and LRS2) to construct general Darwin-Riemann equilib- 
rium models. In Sect. 2.1 we briefly summarize the basic as- 
sumptions and the equilibrium relations are derived in Sect. 2.2. 
More details about the method in general, as well as the applica- 
tions to compact objects, can be found in LRS1 and LRS2 and 
references to equations in these papers are denoted with num- 
bers preceded by "I" and "II", respectively, in the present pa- 
per. Solutions to first order in the deformation will be derived 
for tidal and rotational deformations in Sect. 3.2 and 3.1, respec- 
tively. 

2.1. Model description 

Consider an isolated, self-gravitating fluid system in steady state. 
The system is characterized by conserved global quantities such 
as its total mass M and total angular momentum /. The basic 
idea in our method is to model our self gravitating system by a 
limited number of parameters x\, x>, in such a way that the 
total energy can be written, 



E = E(x u x 2 ,...; M, /,...), 



(1) 



An equilibrium configuration is then determined by extremizing 
the energy according to 



dE 



= 0, j = 1,2, ... 



(2) 



An expression like Eq. (1) can be written down for the total 
energy of a binary system (with components of mass M and M'). 
We adopt a polytropic equation of state between the pressure P 
and the mass density p, 



P = Kp 



1 + 1/n 



(3) 



This defines the polytropic index n and the "entropy" K - both 
are constant within the object and sufficient (with M) to describe 
the mechanical structure of a given object. Under the combined 
effects of centrifugal and tidal forces, the objects (stars or plan- 
ets) achieve nonspherical shapes. We model these shapes as tri- 
axial ellipsoids of principal axes (ai, 02,^3) and (d v a' 1 ,a'^), re- 
spectively. Throughout this paper, unprimed quantities refer to 
the component of mass M while primed quantities refer to the 
component of mass M' . The three directions along which our 
principal axes are measured are, respectively, the line connecting 



the center of mass of the two components, its normal contained 
in the orbital plane and the direction of the orbital angular mo- 
mentum vector. In the simple case of coplanar and synchronous 
rotation, aj is simply the polar radius and ai and 02 are the equa- 
torial radii of the component measured toward its companion and 
in the orthogonal direction, respectively. 

Specifically, we assume that the surfaces of constant den- 
sity within each object can be modeled as self-similar ellipsoids. 
The geometry is then completely specified by the three princi- 
ple axes of the outer surface. Furthermore, we assume that the 
density profile p(m) inside each component, where m is the mass 
interior to an isodensity surface, is identical to that of a spher- 
ical polytrope with the same volume. The velocity field, v, of 
the fluid is modeled as either uniform rotation (corresponding to 
the case of a synchronized binary system), or uniform vorticity, 
V x v, (for nonsynchronized systems). The vorticity vector is 
assumed to be everywhere parallel to the orbital rotation axis. 

For an isolated rotating gaseous sphere, these assumptions 
are satisfied exactly when the fluid is incompressible (polytropic 
index n = 0), in which case the true equilibrium configuration 
is a homogeneous ellipsoid (Chandrasekhar 1969). For a binary 
system, our assumptions are strictly valid in the incompressible 
limit only if we truncate the tidal interaction at the quadrupole 
order. We adopt this quadrupole-order truncation of the interac- 
tion potential in this paper. 

Adding up the orbital separation, r, our set of unknowns 
is (r, 01,02,0?,, a'^a'^a'^) or equivalently (r,p c , A\, A2,p' c , A[, A' 2 ) 
where p c is the central density and A\ = (03 /a i) 2 ^ 3 , A2 = 
(a3/«2) 2 ^ 3 (when no ambiguity exists or otherwise stated, primed 
quantities are defined in the same manner as unprimed ones by 
simply making the transformation x ±± x' throughout the equa- 
tions). The total energy can be written 



E = U + U' + W + W + T + Wi, 



where 



/nP 
— dm = kiKp l J"M 

is the internal energy of component 1 (cf. Eq. (1.3.1)), and 

3 GM 2 



W = -- 



5-n R 



-f = -k 2 GM 5/3 pl /3 f 



(4) 



(5) 



(6) 



is the self-gravitational energy of component 1 (cf. Eq. (1.4.6)) 
with, 



n(n + 1) 

h = — c — -fiie'il, 



5-n 



T 3 (4n\ffi\\ 



1/3 



f(AuA2) = 



A\a\ + A2a 2 + ^3 fl 3 

2(flifl2a3) 2 ^ 3 



A; = Qifl2 fl 3 



Jo fL 



du 



(a 1 + u) 



(7) 



(8) 



(9) 



(10) 



and 



II 2 = (a 2 + u)(a\ + u){a\ + u). 
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£1 and 6\ are the classical variables (dimensionless radius and 
density taken at the surface) used to describe polytropic gaseous 
spheres (see Chandrasekhar 1939) and G is the gravitational con- 
stant. U' and W are similarly defined. The kinetic energy in the 
inertial frame reads 



T = T S + T' + T 



(11) 



where the spin kinetic energy of body 1 (T s ) is given by (cf. 
Eq. (1.5.6)) 

1 2 
T s = -/(A 2 + Q 2 ) - -K n Maia 2 An, (12) 

with Q being the rotational orbital velocity, A is a measure of 
the internal rotation rate in the co-rotating frame, 



5 jf'yg 4 ^ 

3 



(13) 



is a dimensionless coefficient measuring the inertia of the body, 
and 



1 



-K n M{a\ + a\) 



(14) 



is the moment of inertia with respect to the rotation axis. Tables 
giving values of the polytropic constants k\, k 2 , k„, and £i 
as a function of n can be found in LRS 1 and in Chandrasekhar 
(1939). For non synchronous rotation (A + 0), our gaseous body 
is not in the state of solid-body rotation. A rotation rate can 
thus not have the usual meaning. To have a sense of the angular 



velocity, one can take the half of the vorticity (a> 



where 



v is the fluid velocity vector in the inertial frame) as a proxy 
is related to A by 



Q 



9 9 
a \ + a 2 

2a\a 2 



A. 



The orbital kinetic energy Tq is simply 



To = z 



1 MM' 



2M + M' 



2.2 



(15) 



(16) 



Finally, the gravitational interaction energy W\ reads 



GMM' GM , 
Wi = ^(2/1,-4-4) 



GM' 



2r- 



-(2/n - I 22 - I33), 



with 



h = 



Ml - -^KnMa^iy 



(17) 



(18) 



1 This choice has no impact on the result to first order because i) there 
is no cross correlation between tidal and centrifugal distortion at this 
level and ii) the value of u> only plays a role to compute the rotational 
distortion for which solid body rotation is ensured because ci[ = a^. 
In this case, the half vorticity reduces to the usual rotation rate and 
<±> = Q - A. To higher order, this simply highlights the absence of a 
solidly rotating state and the inadequacy of the parametrization by a 
rotation rate in such cases. 



2.2. Equilibrium relations 

We can now derive the set of equilibrium relations 



dE 







M,J,... 



i=l,...,7 



yielding seven algebraic equations for our seven unknowns 
(r,p c , A\, Ai,p' c , A' v AL). The details of the transformations nec- 
essary to express the total energy as a function of the unknowns 
and conserved quantities alone and to be able to carry out the 
differentiation are explained in §2.2.1 of LRS2 and just add tech- 
nical details not needed here. We will thus give only the results. 
Differentiation with respect to r simply yields the modified 
Kepler's law for the orbital mean motion Q 



Q 



G(M + M') 



(1 +A + A'), 



with 



A = 



3 (2/11-/22-/33) 
2 Mr 2 



(19) 



(20) 



Differentiation with respect to the central density p c yields 
the virial relation, 



-U + W + 2T, 



GMM' 



n R 
with the mean radius R = (aiG^c^) 1 / 3 and 



- R m 



2 R 

hi - I33) = r -A. 

3 r 



(21) 



(22) 



Using expressions for U and W, we get the equilibrium mean 
radius 



R = R 



f(A u A 2 ) 1-2 



'\W\ 



5-n\ M' 



M 



gt 



(23) 



where Rq is the radius of the unperturbed spherical polytrope 
given by (Chandrasekhar 1939) 



(n + l)K 
AnG 



(24) 



Finally, the differentiation with respect to Ai and A 2 yields after 
some algebra (cf. Eqs. (1.8.4), (1.8.5) and (1.8.6)): 



2 + 



a 2 AQ A 



2 



a, + a- 



(a\Ai 



qnP-R 

2a 1 AO. A 2 

H 



■ ajA 3 ), 



«2J"R 



2 2 
a 2 + fl 3 



where 



q,,HR 



q n = K„(l - -), 



— {a L 2 A 2 - a^Ai), 



(25) 



(26) 



pn = GM' jr , pr = pn/inGp), and p = M/{%nR ) is the mean 
density of the ellipsoid. 
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3. The shape of gaseous bodies 



3.1. Rotational deformation: Maclauhn spheroids 



In this section, we derive analytical expressions for the defor- 
mation - induced either by centrifugal or tidal potential - of a 
gaseous body. To test the validity of our assumptions, we com- 
pare our predictions with the measured values for the major plan- 
ets of our solar system in Sect. 3.1. Since we do not make any 
assumption about the masses of the two components, these equa- 
tions can be used indifferently to compute the shape of the star 
or of the planet by choosing M = M p and M' = M* when con- 
sidering the planet, and vice versa when focusing on the star. 

In general, the set of equations described in the previous 
section must be solved numerically, but we will study here the 
first order development of these equations at large orbital sepa- 
ration. This approximation corresponds to neglecting terms of 
order OiR^/r 5 ), which is consistent with our truncation of the 
gravitational potential at the quadrupole order and is appropriate 
to address close-in transiting planetary systems. In practice, this 
is done by setting A = A' = and 



fli = R (l + QTi), 



(27) 



in Eqs (19), (23), (25), (26) and their primed equivalent and by 
expanding these equations to first order in a; ((9(Rg/r 3 )). First 
we derive some general formulae by expanding the integrand in 
the definition (10): 



A ; = Rl(l + ct\ + a 2 + a 3 )x 



o 

which yields 



Rq + u 



(a\ + a 2 + a 3 + 2ai) 



du 



2 4 4 , 

At = - + — (ari + a 2 + a 3 ) - -a; + 0{af ). 



To first order, 



and 



f(A u A 2 ) = l+0( ai 2 ), 



A l a ] -A 3 a 3 = — (a x - a 3 )R Q , 
A 2 a\ -A 3 aj = y^<>2 - a 3 )Rl. 



(28) 



(29) 



The principal moment of inertia of the body can also be com- 
puted and reads 



1 



/ =-K n M{a\ + a\) 
2 9 

>-k„MRo(1 +ai +a 2 ) 
2 , 

--k„MRq (spherical case). 



(30) 



The other moments of inertia can be computed by replacing 1 
and 2 by the appropriate indices. The dimensionless moment of 
inertia k„ for different planetary masses, age and stellar irradia- 
tion can be found in tables B.2 and B.5. 



Our set of equations also allows us to compute the effect of the 
centrifugal force alone on a slowly rotationg fluid object. To do 
so, one just has to take the M' -> limit in Eqs (19), (23), (25) 
and (26). Therefore, Q is a free parameter (the rotation rate of 
our body) and there is a degeneracy between Q and A that allows 
us to choose A = 0. 

We introduce the dimensionless angular velocity 

-2 _ <l}lR l nn 

U = (31) 

GM 

as a small parameter of order 0(a\) in all expansions. The vol- 
ume expansion factor can be calculated using 



\W\ 3 GM 2 3 
5-n R 



q n U) 2 . 



(32) 



We get 



\ — a\ + a 2 + a 3 



In 
3-n 



q n co 



(33) 



To the same order of approximation, the two remaining equa- 
tions are given by Eqs. (25) and (26) which yield 



q„(D = -(ai, 2 - a 3 ) 



(34) 



Combinations of Eqs. (33) and (34) give the three figure func- 
tions [cf . Eq. (A 1 2) of LRS2] 



a\ = a 2 — 



15+n 



q n cj 



43 -n 
1 5-3« 



a 3 = ~- 



-q„oj 



2 3-n 

For this configuration, the usual variables are the oblateness 
a\ — a 3 5 



(35) 



/ = 



= a\-a 3 - -q„ti> 
ci\ 4 



(36) 



and the dimensionless quadrupole moment of the gravitational 
field J 2 given by the theory of figures to first order (Zharkov 
et al. 1973) 



j = If-^L = 5g« ~ 2 u 2 Rl 
2 3 J 3 6 GM 



(37) 



For practical purposes, Rq can be computed to first order by us- 
ing Eq. (27) and reads 



«b = J3 



CL\CL 2 Cl 3 

a\ + a 2 + a 3 



\ 27?, 



3 ^eq^pol 



eq ^pol 



(38) 



in our geometry, where R eq and R po \ denote the usual equatorial 
and polar radii. For an incompressible body (n = 0) we retrieve 
the usual solution of the theory of planetary figures / = 5/4 to 2 
(Zharkov & Trubitsyn 1980). 

Attempts have been made to constrain the oblateness and 
thus the rotation period of transiting planets by using the solar 
system planets as test cases (Carter & Winn 2010a,b). Because 
of the wide variety of exoplanets, it is important to have the abil- 
ity to predict the flattening of fluid planets for a wider range of 
parameters than encountered in the solar system. Fig. 1 shows 
the predicted oblateness for various planet masses as a function 
of the rotational period P lot = — . 
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Figure 1. Oblateness given by Eq. (36) as a function of the ro- 
tation period (in days) at 1 Gyr for planets of mass: 0.3Mj (dot- 
ted), 0.5Mj (dashed), 1M, (dash-dotted), 3Mj (long dashed), 
15Mj(solid). The oblateness decreases when the mass of the 
planet increases because massive objects are more compressible 
(see §4), have a more intense self-gravity field and are thus less 
subject to perturbations. 



3.2. Tidal deformation: 
determination of the Love number (k 2 ) 

To compute the shape induced by the tidal force alone, we con- 
sider a non-rotating configuration (to = 0). From Eq. (15), this is 
achieved if A = 2aia 2 £l/(a 2 + (&)■ Then 



T s =-/Q 2 
2 



1 - 



a\+a\, , 



= -m 2 aia 2 = O + CW) 

Thus from Eq. (23) we see that there is no change of volume to 
lowest order, 

fllfl2fl3 



-1=0. 



Since Q 2 = // R (l + p) with p = M/M' and /2 R = O^/r 3 ), only 
the zeroth order must be taken in the left hand side of Eqs. (25) 
and (26), which yields (with help of Eq. (29)) 



a 2 ■ 



15 1*| 
A qn pr^ 



(39) 



Thus [cf. Eq. (A25) of LRS2] 



Ct2 = ai = 



5 lR 3 Q 
2 p r J 
5 IR[ 

—79n~- 

4 p r 



(40) 



As long as the hydrostatic equilibrium holds, this equation can 
be used to compute the shape of the planet and its host star at 
each point of the orbit. We recover the usual dependence of 

the tidal deformation in -j^r, with a factor of order unity, q n , 
which encompasses all the structural properties of the gaseous 
configuration. 



Since we are in the linear approximation with a gravitational 
potential restricted to quadrupolar order, the shape of our body 
can be described with the usual Love number of second order, 
ki (Love 1909, which is twice the apsidal motion constant often 
called £2 in the stellar binary literature). Indeed, once ki (and 
h.2 = 1 + &2 for a body in hydrostatic equilibrium) is known, 
the external potential and the shape that a body will assume in 
response to any perturbing potential can be computed as detailed 
in Appendix A. To derive ki, we compute the quadrupolar term 
of the gravitational potential energy of the system formed by our 
compressible ellipsoid and a point mass, by introducing Eq. (40) 
in the linearized version of Eq. (17), and identify this term to the 
potential energy due to tides given by (Darwin 1908) 



= -k 



GMM'Rq 



(41) 



This yields 



2 1 - f 



-i*H)- 



(42) 



As expected, in the n — limit, we retrieve the Love number of 
an incompressible ideal fluid planet k 2 = 3/2. We can also see 
that k 2 is linked to the square of the dimensionless moment of 
inertia K n . This is because level surfaces are self-similar in our 
model and that the love number encompass both the deformation 
of the body (oc k„) and the gravitational potential created by the 
deformation (oc k„). The Love number for different planetary 
masses, age and stellar irradiation can be found in tables B.2 
and B.5. 

We can see that the value of the Love number tends to de- 
crease with mass above lMj. This is due to the fact that more 
massive objects are more compressible and thus more centrally 
condensed (See Sect. 4). At constant mass, enrichment in heavy 
elements toward the center (possibly in a core) acts to decrease 
the value of k%. In general, redistributing mass from the external 
to the internal layers, which are less sensitive to the disturbing 
potential, decreases the response of the body to an exciting po- 
tential, which translates into a lower k 2 . 

Our model predicts k 2 values in the range 0.3-0.6. As dis- 
cussed by Ragozzine & Wolf (2009) such values of the Love 
number could be inferred by the measurement of the precession 
rate of very Hot Jupiters on eccentric orbits. Such measurements 
could be carried out by Kepler for WASP-12b analogs with an 
eccentricity > 3 x 10~ 4 (most favorable case) or Tres-3 b analogs 
with an eccentricity > 2 x 10~ 3 (for k 2 ss 0.3) and lower ec- 
centricities for higher Love number values. Such measurements 
would indeed be extremely valuable as they would put direct 
constraints on the central enrichment in heavy elements inside 
close Hot Jupiters, like the measurements of the gravitational 
moments of the solar system planets. 

3.3. Synchronized planets 

For values of the tidal dissipation factors inferred for Jupiter 
(Goldreich & Soter 1966; Leconte et al. 2010), the timescale 
of pseudo synchronization of close-in giant planets is less than 
about a million years. The planet is thus in a state of pseudo syn- 
chronization, with a rotation rate given by (Hut 1981; Leconte 
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et al. 2010) 



1 + 



f * 4 + h 6 



(l+3e 2 + |e 4 )(l-e 2 ) 3 / 2 



Q 



(43) 



in the weak friction theory, with e being the eccentricity of the 
orbit. For the simple case of a circular orbit, the spin is thus 
synchronized and, either solving Eqs(19), (23), (25) and (26) 
in the synchronized case (A = 0), or simply adding the results 
of Eqs. (40) and (35) (there is no cross correlation terms to first 
order) yields 



1 



\ +P \Rl 

.3 



U2 = 



Qf3 = 



-3*. 



-3* 

















^0 






r3 





l+p 
l + p 

2-P 
l+p 

5 + 2p 

T+~p~ 



In 
In 



3-n 
In 



3-n 



and 



R 3 



1 +q,, 



In 
3-n 



i + p\K 

.3 



(44) 



(45) 



3.4. Model Validation 

There are two major assumptions in the present calculations: 

- The absence of a central core. The aim of such an approx- 
imation is to avoid to introduce any free parameter in the 
model. In any case, the core mass and the global enrich- 
ment in giant extrasolar planets are yet weakly constrained 
(Guillot 2005; Leconte et al. 2009). We will show that this 
approximation introduces an uncertainty of w 10% on the 
derived shape. 

- The polytropic assumption. This allows us to derive a com- 
pletely analytical model. Comparison with a more detailed 
numerical integration (See Appendix A) shows that the de- 
viation between the results of the two models (analytical vs 
numerical) is smaller than the uncertainty due to the no-core 
approximation. 

Since the oblateness (/), J2, hi, the mean radius and rota- 
tion rate are known for the major planets of our solar system, we 
can test our theory on these objects. The details of the calcu- 
lation of the chosen polytropic indexes are presented in Sect. 4. 
The results are summarized in table 1, which shows the actual 
values of the relevant parameters for the two major planets taken 
from Guillot (2005), the values of the oblateness and of J2 calcu- 
lated with our model and, for comparison, with the assumption 
of an incompressible body (« = 0). We see that, whereas the 
values of / derived from the incompessible model differ from 
the true values by almost a factor of 2, our polytropic model 
predicts the /-values to within 12%. Note that higher-order 
terms (of order at 2 ~ 10%) are not completely negligible for 
rapidly rotating bodies such as Jupiter and Saturn (see Zharkov 
et al. 1973; Chabrier et al. 1992). The polytropic model, how- 
ever, yields the /2-values that differ from the measured values 
by 30% (for Jupiter) and 59% (for Saturn). These discrepan- 
cies are mostly due to the large metal enrichment in these planet 
interiors, probably with the presence of a large dense core as de- 
tailed three layers models can reproduce exactly the measured 



moments (Chabrier et al. 1992). Note that this no-core approx- 
imation has less relative impact on the distortion of the shape 
predicted by the model than on the gravitational moments be- 
cause these effects scale as hi — 1 + ki (See Appendix A) and 
lc2 respectively, Ic2 being < 0.6 in the situations of interest. Such 
discontinuities in the density profile (and its derivatives) could 
be addressed more precisely with two different polytropes, but 
this would add extra free parameters and would not serve the 
very purpose of the present paper. 

Table 1. Comparison between the measured oblateness and data 
of the gravity fields and the values obtained with our polytropic 
model (polytrope) and with a model with n = (incompress- 
ible). 





Jupiter 


Saturn 


M p [10 2b kg] 


18.986112(15) 


5.684640(30) 


tf eq [10 7 m] 


7.1492(4) 


6.0268(4) 


fi P oi[ lo7 m] 


6.6854(10) 


5.4364(10) 


R [10 7 m] 


6.9894 


5.8198 


P mt [10 4 s] 


3.57297(41) 


3.83577(47) 


or x 10 2 


8.332 


13.940 


h 


0.49 


0.32 


n 


0.936 


0.748 


q n 


0.547 


0.623 


fx 10 2 


6.487(8) 


9.796(9) 


fx 10 2 (polytrope) 


5.701 


10.849 


/ x 10 2 (incompressible) 


10.416 


17.425 


J 2 x 10 2 


1.4697(1) 


1.6332(10) 


J2 x 10 2 (polytrope) 


1.023 


2.586 


Ji x 10 2 (incompressible) 


4.166 


6.970 



The numbers in parentheses are the 
uncertainty in the last digits of the value. 



While we decided to use a polytropic assumption to infer 
a fully analytical theory, the figures of a body in hydrostatic 
equilibrium can be derived without this assumption. As shown 
in Sterne (1939) (See also Zharkov et al. 1973; Zharkov & 
Trubitsyn 1980; Chabrier et al. 1992 for more detailed appli- 
cations to the giant planets case) and outlined in Appendix A, 
this theory, however, requires a numerical integration even to 
first order. For an ideal n = 1 polytropic sphere, Eq. (42) agrees 
with the numerical results of Sterne (1939) with less that 1% 
error. To compare these methods in our context, we derive the 
values of Ic2 using our analytical model (Eq. (42)), and by numer- 
ical integration of Eqs. (A. 6) and (A. 11) for our best representa- 
tive models of Jupiter and Saturn in our grid (Although without 
cores). For Jupiter, our Eq. (42) gives Ic2 = 0.55, the numerical 
integration gives &2 - 0.57 to be compared with the measured 
value of &2,j = 0.49. For Saturn, these values are 0.60, 0.66 
and ^2,s = 0.32 respectively. Both models predict Ic2 values than 
are higher than the measured ones. A direct consequence of the 
presence of heavy elements inside our giant planets. Comparing 
the models, our Eq. (42) yields slightly smaller £2 values than 
the numerical integration which tends to mimic a central over- 
density (See Sect. 3.2). As discussed above a more precise mod- 
eling requires the addition of central enrichment in heavy ele- 
ments whose mass fraction would be a free parameter. Without 
better knowledge of the internal composition of giant exoplan- 
ets, we think that the two methods yield similar results up to the 
sought level of accuracy. For sake of completeness, the values of 
k.2 computed with both methods are presented in tables B.2 and 
B.5. 
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4. Polytropic index in gaseous irradiated planets 

To readily use the results of Sect. 3, one only needs to have a 
proper value for the polytropic index n to be used. In this sec- 
tion, we derive realistic polytropic indices from numerical mod- 
els of gaseous irradiated planets. All the other polytropic func- 
tions (k„, q„, ...) can be derived by integrating the Lane-Emden 
equation and are tabulated in Chandrasekhar (1939) and LRS1. 
They are given for different planetary masses, age and stellar 
irradiation in tables B.2 and B.5. We focus on the polytropic 
index in the planet because, in the context of transiting exoplan- 
ets, both the stellar rotation and the stellar tides have a negligible 
impact on the transit depth, as will be discussed in Sect. 5. The 
main physics inputs (equations of state, internal composition, 
irradiated atmosphere models, boundary conditions) used in the 
present calculations have been described in detail in previous pa- 
pers devoted to the evolution of extrasolar giant planets (Baraffe 
et al. 2003; Chabrier et al. 2004; Leconte et al. 2009), and wiU 
not be repeated here. 

We computed a grid of evolution models of gaseous 
giant planets with solar composition for various masses 
M p e [0.35Mj, 20Mj] and incoming stellar flux 
F+ e [0,4.18 x 10 6 W.irT 2 ]. Low irradiation model 
can be used to infer the oblateness of long period rotating 
planets (such as Jupiter), while strongly irradiated models 
can be used to infer the shape and its impact on the transit of 
close-in planets. For the non irradiated case, the grid extends 
to 75Mj. As the effect of irradiation on the internal structure 
decreases with the effective temperature of the object, these 
models computed with non irradiated boundary conditions 
should give a fair description of massive brown dwarfs in the 
range of irradiation considered. The pressure-density profile 
of each model is then fitted by a polytropic equation of state 
(Eq. (3)) at each time step and an example of the result of such 
a fit is shown in Fig. 2. Note that the disagreement between the 
actual P - p profile and the polytrope in the lower left area of 
Fig. 2 is both expected and needed: This low-density region (the 
first 5% in mass below the atmospheric boundary surface) has a 
different effective polytropic index than the planetary interior. 
In order to capture the bulk mechnical property of the planet, 
we weight each shell in the internal structure profile by its 
mass during the fitting procedure. This provides us with a grid 
tabulating the polytropic index of the planet, n p = n p (M p , ,F*, f), 
and its spherical equilibrium radius, Rq p = Rq p (M p , F+,t), 
where t is the age of the object. These functions, along with 
other quantities ( T e ff, ...), are tabulated in tables B.2 and B.5 2 . 
Figures 3 and 4 show the variation of n with the mass for 
different ages with F+ — and 4.18 x 10 6 W.irT 2 , respectively. 

As shown in Fig. 3, in the non-irradiated case, we re- 
cover qualitatively the results of Chabrier et al. (2009): ex- 
cept for the early stages of the evolution, the (dimensionless) 
isothermal compressibility of the hydrogen/helium mixture is 
a monotonically increasing function of the polytropic index, 

X = f^f = if^> an d thus of the mass of the object. In the 
high mass regime, n slowly increases as the relative importance 
of ionic Coulomb effects compared with the degenerate electron 
pressure decreases, and approaches the n = 3/2 limit, the ex- 
pected value for a fully degenerate electron gas, when M p ap- 
proaches the hydrogen burning minimum mass (« 70Mj) as can 
be seen in table B.4. In the low mass regime, the compressibil- 



2 Electronic versions of the model grids are avail- 
able at http://perso.eiis-lyon.fr/jeremy.leconte/JLSite/JLsite 
/Exoplanets _Simulations.html 
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P (kg/ m ) 
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Figure 2. The internal pressure-density profile of an irradiated 
1.8Mj planet (solid line). The dashed line represents the best-fit 
polytropic equation of state. The pressure-density range covered 
in the inner part of the body (95% in mass) is represented by 
the thicker part of the solid curve, which is well modeled by a 
polytropic EOS. As the thin part of the P — p curve represents 
only 5% in mass of the body it is disregarded by the fit. 
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Figure 3. Polytropic index for non-irradiated planets as a func- 
tion of the planet's mass M p at 100 Myr (Dotted), 1 Gyr (Dashed) 
and 5 Gyr (Solid). The shaded area represents the uncertainty on 
the polytropic index for the 5 Gyr case (see text). 



ity decreases with the mass because the repulsive Coulomb po- 
tential between the ions, and thus the ionic electrostatic energy 
becomes dominant. Ultimately, electrostatic effects dominate, 
leading eventually to;f » n « for solid, terrestrial planets. 

A new feature highlighted by the present calculations is the 
non-monotonic behavior occurring between 1 - 3 Mj at early 
ages. This occurs when the central regions of the planet, of pres- 
sure P c and temperature T c , previously in the atomic/molecular 
regime, become pressure-ionized, above 1-3 Mbar and 5000- 
10 000K (Saumon et al. 1995; Chabrier et al. 1992; Saumon et al. 
1992), and the electrons become degenerate. An effect more 
consequential for the lowest mass objects, whose interiors en- 
compass a larger molecular region. This stems from the fact that 
(Chandrasekhar 1939) 



P c > 



GMl 
8tt^' 



and 



GM p 



2-3 Mbar * P 



ionization ■ 



(46) 
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Older (with smaller R p ) and more massive (M p > 2Mj) objects 
have P c > 10Pi on i zat i O n and the ionization extends all the way 
up to the outermost layers of the gaseous envelope, which then 
contains too small a mass fraction of molecular hydrogen to sig- 
nificantly affect the value of the polytropic index. This contrasts 
with younger objects around 1-3 Mj, whose external molecular 
hydrogen envelope contains a significant fraction of the planet's 
mass, leading to a larger value of the polytropic index, as molec- 
ular hydrogen is more compressible than ionized hydrogen (see 
e.g. Fig. 21 of Saumon et al. 1995). Once again, for these latter 
objects, the interior structure would be better described by us- 
ing two different polytropes, but such a significant complication 
of the calculations is not needed at the presently sought level of 
accuracy. 

As seen on Fig. 4, a strong irradiation enhances the afore- 
mentioned feature: the evolution is delayed because the irradi- 
ated atmosphere impedes the release of the internal gravothermal 
energy. This yields a slower contraction, thus a lower central 
pressure (and lower central temperature) for a longer period so 
that the object enters the ionization regime at a later epoch. The 
bump at the high mass end of the 100 Myr isochrone is due to 
deuterium burning which also occurs later for a given mass, be- 
cause of the cooler central temperature (see above). At 100 Myr, 
the 20 Mj has already burned a significant amount of its deu- 
terium content and starts contracting again, whereas lower mass 
planets are still burning some deuterium supply, leading to a less 
compact and thus less ionized structure. This leads to the non- 
monotonic behaviour on the high-mass part of the n-M diagram 
at 100 Myr, which reflects a similar behaviour in the mass-radius 
relationship. 



and the central density, using 

Since our numerical simulations provide both the ra- 
dius, /?o,p(Af p ,F*,f), and the central density of the object, 
Pc,p(-^p>-P\t> 0> we can invert Eqs. (24) and (47) to compute K v 
and n p . This new determination of the polytropic index is com- 
pared with the previous one, obtained by fitting the P—p profile, 
in Figs 3 and 4 for the 5 Gyr case: the new n p value corresponds 
to the upper envelope of the shaded area. Fig 3 shows that the 
two approaches yield very similar results in the non-irradiated 
case. For the irradiated case, the average uncertainty on our de- 
termination of « p lies between about 5 % and 15 % for the low 
mass planets. 

5. Implications for transit measurements 

When limb darkening is ignored, the depth of a transit is given by 
the ratio of the planetary and stellar projected areas. When both 
bodies are spherical, this simply reduces to (5L*/L+ oc (R p /R+) 2 . 
For close-in planet-star systems, however, both tidal and rota- 
tional deformations yield a departure from sphericity, so that 
what is measured is no longer the mean radius but an effective 
"transit radius" defined such that the cross section of the planet 
is equal to 7iR 2 and similarly for the star. Thus the transit depth 
6 reads 




10 

M P (M j up ) 



Figure 4. Polytropic index for strongly-irradiated planets as 
a function of M p at 100 Myr (dotted), 1 Gyr (dashed) and 5 
Gyr (solid). As the irradiated atmosphere impedes the radia- 
tive cooling of the objects, it retards its contraction. Therefore, 
the non-monotonous behavior observed at the early ages in the 
non-irradiated case (Fig. 3) is enhanced, even at a later epoch. 
The bump at the high mass end of the 100 Myr curve is caused 
by deuterium burning (see text). The shaded area represents the 
uncertainty on the polytropic index for the 5 Gyr case (see text). 



To evaluate the uncertainty in our determination of the poly- 
tropic index, we use an alternative method to derive n. As shown 
by Chandrasekhar (1939), the knowledge of M, K and n is suffi- 
cient to infer the radius of the polytrope, with the help of Eq. (24) 



5. 1 . Impact on transit depth 

In general, the projected area of an ellipsoid can be computed for 
any orientation and then at each point of the orbit, as explained 
in Appendix B. Figure 5 shows the projected area of the planet 
(71R 2 ) as a function of its anomaly (<p) and inclination (i) 



R 



tr,p 



tJci 2 sin 2 i (t 



a 2 sin 2 1 



+ a 2 cos 2 1 



+ a\a\ cos 2 i, 



(49) 



normalized to the spherical case (7iR ( 2 p ). 

When the planet is seen from its "side" (<f>/n = 0.5), the 
observer sees a bigger planet because the rotation of the latter 
on itself tends to increase its volume, as has been mentioned by 
Li et al. (2010) for WASP-12 b. The possibility to measure these 
effects from the light curve is discussed in Ragozzine & Wolf 
(2009) and Carter & Winn (2010a). 

For the simple case of an edge-on orbit at mid transit (<p = 0, 
i = 90°), since the observer, the planet and the star are aligned 
with the long axis of the tidally deformed ellipsoid 3 ' 4 , R tr p = 
V fl 2,pfl3, P and R tu + = ya 2; *a3,*- Therefore, 



S= W* . (1+77)) 



(50) 



where Ro p and Rq ^ are the respective radii the planet and the star 
would have in spherical equilibrium and rj is by definition the 
variation of the transit depth induced by the ellipsoidal shape 



3 This is still verified to first order in ( 
terms appear. 



and i — | as only second order 



4 In the following, the variables have the same meaning as in Sect. 2 
and 3 with p indices when referring to the planet and ★ to the star 
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Figure 5. Normalized projected area of the planet as a function 
of its anomaly ((/>) for inclinations of the orbit going from i = 90° 
(lowest curve) to i - 0° (highest curve) by steps of 10° (5° for 
bottom panel) for a WASP-12b analog on a circular orbit. Top: 
For the full orbit. Bottom: zoom on the (primary or secondary) 
transit. The ordinates of the dotted, solid and dashed horizon- 
tal lines are respectively a\a 2 IR\ (face-on orbit), a\a^jR\ and 
a 2 a 3 /Rl. 



of the components relative to the transit depth in the spherical 
case. To first order in the deformation, this is given by 



T] = ar 2 ,p + <*3,p - 012,* ~ Ofy*. 



(51) 



The choice of the expression to be taken for the or} depends on 
the physical context. In the general case, one can use a linear 
combination of Eqs. (35) and (40) and get a general expression 
which depends on r, u p and <y+. However, most of the planet 
hosting stars have a low rotation rate compared to the orbital 
mean motion. This entails that the rotational deformation is neg- 
ligible compared to the tidal one and can generally be neglected. 
As mentioned in Sect. 3.3, hot Jupiters should be pseudo syn- 
chronized early in their evolution. Therefore, we will assume 
such an approximation in our calculations in order not to intro- 
duce any other free parameter. The impact of the rotation alone 
is described Sect. 3.1. Under such an approximation, 



1 



l+p 



0,* 



0,p 



511 + P 
4\l+p 



4« p 
3 - n r 



2 r J 



(52) 



where the parameter p now denotes the mass ratio M p /M*, and 
q p and q+ are equal to q„ for n = n p and n — n+, respectively. 



The first line in the above equation represents the contribution 
of the planet, which is always negative (for reasonable values of 
n). Our line of sight follows the long axis of the tidal bulge and 
we see the minimal cross section of the ellipsoid. 

The contribution of the star is positive and, in most cases, 
negligible compared the planet's contribution because 



R 3 n 1 + p 

0,p 1 



«: 1, 



for a typical system (1CT 3 for a Jupiter-Sun like system). As a 
consequence, the results presented hereafter do not depend on 
q+ as long as realistic values of n* e [1.5, 3] are taken. 

Figure 6 portrays the relative transit depth variation com- 
puted with Eq. (52) for several planet masses as a function of 
the orbital distance, for a Sun-like parent star. While all the 
curves are calculated at an age of 1 Gyr, they do not change 
much for older ages because both the radii and the polytropic 
indices remain nearly unchanged after 1 Gyr (see Fig. 4). Given 
the accuracy of the radius determination achieved by the latest 
observations (1 to 10%), the transit depth variation is significant 
for Saturn mass objects (M p w Mj/3) closer than 0.04 AU and 
Jupiter mass objects closer than 0.02-0.03 AU. Because we de- 
rived the equations to first order, the value of tj derived from our 
model should be taken with caution when rj > 0.1 - 0.3 (and are 
clearly not meaningful for 77 > 1). In this regime, correspond- 
ing to the upper left region of Fig. 6, one should use the theory 
of planetary figures to higher order, but then numerical calcu- 
lations become necessary, loosing the advantage of our simple 
analytical expressions. Figure 6 also displays the transit depth 
variation computed for the most distorted known transiting exo- 
planets, with the observationally measured parameters. The er- 
ror bars reflect the uncertainties in the model and in the measured 
data. 



5.2. Which radius? 

Before going further, it is important to summarize the differences 
between the various radii that we have defined above. Note 
that, in the literature, the term "radius" is used loosely, even 
for non-spherical objects. Importantly enough, this can lead to 
discrepant normalizations throughout different studies and pub- 
lished values of transit radius measurements when, for example, 
radii are shown in units of Jupiter radii (Rj) without precisely 
defining the latter. 

One can define a\, a 2 and 03 as the distances between the 
center and some isobar surface along the three principal axes of 
inertia. For any distorted object, we can define the mean radius 
(R) as the radius of the sphere that would enclose the same vol- 
ume as the described surface. In our case of a general ellipsoid, 
we have R = (a\a 2 a-s) 1 ^ . If axial symmetry holds (e.g. for a 
rotating fluid body), we have a\ — a 2 = R eq , defining the equa- 
torial radius, and 03 = R po \ the polar radius. Finally, Rq is the 
radius of the spherical shape that the fluid body would assume if 
it was isolated and at rest in an inertial frame (the limiting case 
for which all the mentioned radii would be equal). This latter 
is the radius computed in usual ID numerical evolution calcula- 
tions. Note that in general R + Rq because the centrifugal force 
has a net outward component that increases the volume of the 
object, as can be seen from Eq. (33). 

One must be aware that only a\, a 2 and 03 (reducing to 
R sq and R po \ for solar system gaseous bodies) can be measured 
directly and are not model dependent. This is why we de- 
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Figure 6. Relative transit depth vari- 
ation rj computed with Eq. (52) as a 
function of the semimajor axis at 1 Gyr 
for planets of mass: 0.3Mj (dotted), 
0.5Mj (dashed), IM, (dash-dotted), 
3M ; (long dashed), 15Mj(solid). 
The shaded area shows the zone 
where higher order terms become 
non-negligible. The decrease of the 
transit depth due to tidal interactions 
is smaller when the mass of the planet 
increases because massive objects 
are more compressible (see §4) and 
thus less subject to non-spherical 
deformations. 



fine Rj as the equatorial radius of Jupiter at the 1 bar level 
(Ri = Reqj = 7. 1492x 10 7 m, Guillot 2005 and reference therein). 

Unfortunately, transit measurements only give access to the 
projected opaque cross section of the planet (= ttR^ v ) defining 
a "transit radius" which depends on the shape of the planet, its 
orientation during the observation and the wavelength used. To 
convert this transit radius inferred from the observations (<R t r,p) to 
the spherical radius (/?o,p) - that can be compared to ID numer- 
ical models - one must eliminate 6 from Eqs. (48) and (50). As 
shown above, the stellar impact on r\ is negligible compared to 
the planet's contribution (^ t r,* ~ ^o,*)- Then, using the first term 
in Eq. (52) and expanding the expression giving the definition of 
/?tr,p, one gets 

R0,p * fltr.p ■ (l - ^(«2,p + OS A * fltr.p • (l - |) . (53) 

For the most distorted known planets, the relative variation be- 
tween the transit radius and the equilibrium radius 

AR = (/? 0) p - R tr , P )/R tr , P * -17/2 

is positive and amounts to 3.00% for WASP-12b, 2.72% for 
WASP-19b, 1.21% for WASP-4b, 1.20% for CoRot-1 b, 0.89% 
and OGLE-TR-56 b. 5 

Note that because the mean density scales as R^ 3 , the in- 
crease in radius implies a decrease in the mean density inferred 
which is about three times larger (i.e. m 9% for WASP-12b). 
This is of particular importance when one wants to constrain the 
internal composition or enrichment of giant planets from transit 
measurements. 



5 Of course, since r\ oc \ -^r, Eq. (53) is an implicit equation on 
/?o, p . To obtain R p to the sought accuracy, a perturbative development 
in powers of n tr = rj(R 0p = R tlp ) can be obtained using recursively 
Eq.(53) ' 

^,1-^!-! + ?£_!ji +O0jt) . (54) 

Ktr.p I I 4 I 

However, terms of order rj^ are of the same order than the second order 
corrections to the shape that we have neglected throughout. 



6. Conclusion 

Because of the large variety of exoplanetary systems presently 
discovered, with many more expected in the near future, and 
the increasing accuracy of the observations, it is important to 
take into account the corrections arising from the non-spherical 
deformation of the planet or the star, due to rotational and/or 
tidal forces, as such a deformation yields a decrease of the tran- 
sit depth. In order to do so, it is extremely useful to be able to 
compute analytically the shape of planets and stars in any con- 
figuration from the knowledge of only their mass, orbital separa- 
tion and one single parameter describing their internal structure, 
namely the polytropic index, n. Such formulae are derived in 
Sect. 3, and can be easily used to determine the impact of the 
shape of the planet on its phase curve and on the shape of the 
transit light curve itself (Carter & Winn 2010a). They can also 
be used to model ellipsoidal variations of the stellar flux that are 
now detected in the CoRoT and Kepler light curves (Welsh et al. 
2010). These formulae also give good approximations for vari- 
ous parameters describing the mass redistribution in the body's 
interior and the response to a perturbing gravitational field, i.e. 
the moment of inertia, /, and the Love number of second degree 

Another major implication of the present work is to show 
that departure from sphericity of the transiting planets produces 
a bias in the determination of the radius. For the closest planets 
detected so far (< 0.05 AU), the effect on the transit depth is of 
the order of 1 to 10% (see Fig. 4), by no means a negligible ef- 
fect. The equilibrium radius of these strongly distorted objects 
can thus be larger than the measured radius, inferred from the 
area of the (smaller) cross section presented to the observer by 
the planet during the transit. The analytical formulae derived 
in the present paper, and the characteristic polytropic index val- 
ues derived for various gaseous planet masses and ages, allow 
to easily take such a correction into account. Interestingly, since 
this equilibrium radius is the one computed with the ID structure 
models available in the literature, the bias reported here still en- 
hances the magnitude of the puzzling radius anomaly (see Fig. 6 
of Leconte et al. (2010)) exhibited by the so-called bloated plan- 
ets. 
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The cos (sin) corresponds to positive (negative) values of m and 
Pj" are the usual associated Legendre polynomials. The refer- 
ence axis defining 6 and if/ may change from one term to the 
other. For example, the rotation axis is best suited to treat ro- 
tational distortion and the line connecting the center of mass of 
each body is better to describe the tidal distortion. It is shown by 
Sterne (1939) that to first order, the shape of the distorted level 
surface of mean radius s (see Zharkov & Trubitsyn 1980 for a 
detailed definition) takes the form 

r(s, 0,iff) = s(l + s|» Y™(6, if/)) , (A3) 

where r is the distance between the center and the level surface 
as a function of 9 and i// and s™(s) a figure function yet to be 
calculated. Sterne (1939) shown that, ignoring terms of order 
s™ X s° , s™ verifies the following differential equation 

i2 



d 2 s™ 1(1+1) 6 p(s) ( *r 



ds 2 



s 2 p(s) \ ds 



with 



* 3 Jo 



f p(s / )s' 2 ds / . 
Jo 



(A.4) 



(A.5) 



Using the variable 771(5) 



s!"(.r) ds 



this rewrites 



+ $ " W " 1(1 + l) + 6 ^ + D = °- ( A - 6 > 
ds p(s) 

Then, J]i(R) (R being the external mean radius of the object) can 
be obtained by numerical integration (with 771(0) = 1-2) and 
the shape and external potential (I/™) of the body are given by, 
respectively 



s m (j>\ - 21 + 1 — c m R 

Sl() " 1 + 7]\(R) MG 



1+1 



and 



1 + m (K) r 1 



(A.7) 



(A.8) 



Appendix A: Theory of planetary figures: 
numerical methods 

Here, we briefly outline the method described in Sterne (1939) to 
compute numerically the response of a body in hydrostatic equi- 
librium to a perturbing potential 6 and derive additional formulae. 
To lowest order (which is consistent with the order of approxi- 
mation used throughout the present paper) the body response is 
linear and the total deformation is the sum of the response to 
each term of the decomposition of the perturbing potential. Let 
us consider a term of the decomposition of the potential of the 
form 



vP(r,0,iW 



where the T, m are tesseral harmonics defined by 



cos |n# p H (cos9l 
sin \m\tfi 1 



(A.l) 



(A.2) 



6 Here we take the convention that the force acting on a particule 
of mass M due to a potential V is F = -MgradV. This yields some 
difference of signs with Sterne (1939). 



In order to compare this numerical model with others, we 
can compute several observable quantities. By definition, the 
potential Love number of degree 2 (k 2 ) is given by 



U° 2 (R, 9, iff) = k 2 V%(R, 0, iff), 



which yields 



k 2 



3 - m(R) 



2 + m (R) 

The level Love number (h 2 ) is given by 

V 2 \R, 6, if/) 



s° 2 (K)Y°(e,ify) = -h 2 



gR 



where g is the surface gravity acceleration, and 

■■ ki + 1, 



2 + m (R) 

as expected for a body in hydrostatic equilibrium. 



(A.9) 



(A. 10) 



(All) 



(A. 12) 
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A.1. Axi-symmetric Case 

Thus the first gravitational moment {J 2 ) defined by 



is given by 



d3 

Ji =k 2 cf 

GM 1 



(A.13) 



(A. 14) 



No distinction is made between R and R eq (the equatorial radius) 
when comparing Eqs. (A. 8) and (A.13) because this would only 
add higher order corrections to J 2 which is already a first order 
quantity. 

For the rotational distortion of the body whose angular ve- 
locity is lo, c9 = u> 2 /3 and 



k 2 oj 2 R 3 
7 GM 



h 



(A.15) 



If one is concerned with the external shape, the oblateness (/, 
see Eq. (36)) of a rotating body is given by 



3 

h 2 co 2 R 3 k 2 + 1 co 2 R 3 



f = -~*l 



2 GM 



2 GM 



(A. 16) 



By extension, one can define J 2 for a tidal perturbation by a 
secondary of mass M' at a distance r', = leading to 



h 



M' (R 



M 



(?) 



(A. 17) 



but the reference axis is the line connecting the two center of 
mass and not the rotational axis. 



A.2. Triaxial Case 

While it is tempting to add Eqs. (A.15) and (A. 17) to obtain the 
total J 2 of a body in a close binary, we must remember that the 
tidal and rotational deformations do not have the same axis of 
symmetry in general. Taking 8 as the colatitude and iff as the 
longitude of the body considered, the external gravitational field 
of the latter reads 

GM 

U(r,9,i//) = X 

r 

00 1 , x i 



1=0 m=0 



Yj Yj I ~T I (C i m cos m( A + 5 ™ sin mtA)P|" (cos ff). 



(A.18) 



The quadrupole moment in the linear approximation, is given at 
the surface by 



Y(C™ cosmt/r + S™ sinmt/0P™(cos6») 

= k 2 R 2 eq (y P 2 \cos 9) - ^P 2 (co S 0')) , (A. 19) 

where ff is the angle between the current point and the line con- 
necting the two center of mass. For the coplanar case where the 



tides raising object orbits in the equatorial plane of the distorted 
body, cos ff = sin 9 cos iff and thus 



n 3 , 3sin 2 (6>) 
P°(cos ff) = - sin 2 (#) cos(2i/r) + 



1 



1 



= - cos(2t/0P^(cos ff) - -P^cos ff). 



Thus 



and 



-a 



1 oj 2 R 3 



eq lM'/R 
+ 



3 GM 2 M 



(*)" 



Ci = -k 2 



1 M' {R\ 3 



-4 M \r 



(A.20) 



(A.21) 



(A.22) 



All the other moments are equal to 0. Similar decompositions 
can be used to infer the precise shape of the surface from a sum 
of perturbing fields. This gives 



r(R, 9, t/r) 
R 



= 1 



1M'/R\ 3 1oj 2 R 3 



(?) 

1 M' IR\ 3 



h 2 1 - — I - I + - "—^ I P?(cos ff) 
\2 M\r') 3 GM ' ; 



+ Kl (I'jf (?) ) C ° S(2tA) P 2 (cos 0) ' (A ' 23) 

which directly translates into a\, a 2 and 03 (once R is known) by 
setting (9, ip) equal to (tt/2, 0), (n/2, n/2) and (0, 0), respectively. 
Translating this into au o: 2 and aj is a little more complicated 
because one needs to account for the fact that R > Rq due to 
the centrifugal potential. This can be taken into account either 
numerically - by including the centrifugal force when solving the 
hydrostatic equilibrium - or analytically using Eq. (32) or (45). 



Appendix B: Projected area of a triaxial ellipsoid 

B.1. General case 

Let us define two coordinate systems. The first one (x', y', z') is 
defined by the three main axes of the ellipsoid. In this frame, the 
equation of the surface of the ellipsoid is 



T 



= 1. 



(B.l) 



To compute the projected area of this ellipsoid as it will be seen 
by the observer, it is easier to put ourselves in another coordinate 
system defined by the line connecting the center of mass of the 
system and the observer (toward the observer; x), the projection 
of the orbital angular momentum on the sky plane (z) and a third 
axis in the sky plane chosen so that (x, y, z) follows the right- 
hand vector sense. The current position vector (r = (x, y, z)) 
expressed in this frame is thus related to the one expressed in the 
first coordinate system by a rotation matrix % such as 



K-r, 



(B.2) 



With 9? tr 9? = 1. The equation of the ellipsoid in the new system 
thus writes 



^(r)Er' r f 



3 ^ 
0^0 

at / 



ftr = r tr yir = 1. 



(B.3) 
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The exact value of the matrix J[ will depend on the rotation 
needed and on the angles chosen to represent it. This can be 
worked out in each specific case. To keep some generality, we 
will take J[ of the form 



M = 



a d f\ 
d b e 
fee, 



(B.4) 



The symmetry is ensured by the fact that both of our coordi- 
nate systems are orthonormal. The equation of the contour of 
the projected shadow is given by the fact that the normal to the 
ellipsoid is normal to the line of sight (x) there. This assumes a 
completely opaque body below the isobar chosen to be the sur- 
face. The complete calculation of the level at which optical rays 
that are grazing, at the terminator, have an optical depth close to 
unity (Hubbard et al. 2001; Burrows et al. 2003; Guillot 2010) - 
in the present geometry - should give rise to subtle effects but of 
smaller importance. This reads 



= grader)]" 
= 2r tr ^x. 



(B.5) 



This shows that these points are located on a plane whose equa- 
tion is (since a + 0) 



x = --(dy + fz). 
a 



(B.6) 



Substituting x in Eq. (B.3) by Eq. (B.6) we see that the cross sec- 
tion is an ellipse following the equation 



(y,z) 



b- d - e-^ 



df 



J 



y z \s(y.:.)>£ 



1. 



(B.7) 



It is thus possible to find the rotation in the sky plane needed 
to reduce the ellipse and find its principal axes (pi, pj). If only 
the cross section {np\p2) is needed, we can use the fact that the 
determinant of a matrix is independent of the coordinate system 
so that 



The matrix can be computed thanks to Eq. (B.3) giving 
the a, b, / coefficients and thus Det(S). This gives the project 
area of the planet or the star at any given point of the orbit 



np 



ip 2 = n sin 2 /( 



a 2 sin 2 i (a 2 sin 2 <p + cos ; 



4>) 



+ a\a\ cos 2 i. 



(B.12) 



as shown on Fig. 5. Note however, that this formula is more 
general. It is also the cross section that would be seen by an 
observer located in the direction 



r obs - 



sin i cos (f> } 
sin i sin <; 
cos i 



(B.13) 



in the reference frame defined by the three main axes of the el- 
lipsoid (x, y, z). 



npipi = 



VDet(S) 



with 



Det(S) = be - e 2 



bf 
a 



cd z 
a 



+ 2 



def 
a 



(B.8) 



(B.9) 



In the case of an edge-on orbit at mid transit, no rotation is 



needed, % is the identity and thus a 
and d — e — f = 0. We retrieve 



l/af, b = l/aj, c 



1/a 2 



y[b~c 



(B.10) 



B.2. Coplanarcase 



If the planet equator and the orbital plane are coplanar, the unit 
vectors of first coordinate system defined above coincides with 
the unit vectors defined by the line connecting the two center of 
mass (from the secondary to the object under consideration; x'), 
its normal in the orbital plane (in the direction of motion; y') and 
the rotation axis of the body (z')- If i is the inclination of the orbit 
with respect to the sky plane and (f> the true anomaly defined to 
be at mid transit, the rotation matrix defined by Eq. (B.2) reads 



n = 



sin(/)cos(0) sin(0) cos(/)cos(0) , v 
- sin(/) sin(0) cos(0) - cos(/) sin(0) 
- cos(0 sin(/) 



(B.ll) 
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Table B.l. Summary of variables 



Symbol Definition 



Equation 



E 
U 
W 
T 
T s 
T u 
Wi 
n, K 
G 
M 
R 
P 

P 
a\ 
ffj 
r 

h 
h 

f 
a, 
/ 

Kn 
Q. 

A 

(-0 
?! 
*i 
h 

A 

gi 
R 
q« 

Mr 
Mr 
P 

/ 
h 

6L< 
8 

R u 

n 



total energy 

internal energy 

self gravitational energy 

total kinetic energy 

rotational kinetic energy 

orbital kinetic energy 

gravitational interaction energy 

polytropic index and temperature 

gravitational constant 

Mass of primary 

Mean radius of the primary 

Pressure 

density 

principal axes of the ellipsoid 
first order correction to a\ 
orbital distance 
asymmetry factor = (a^/ai) 213 
central density 

scaling constant of the internal energy 
scaling constant of the gravitational energy 
ellipsoidal correction to the gravitational energy 

principal moment of inertia 

Dimensionless moment of inertia 

orbital mean motion 

internal angular velocity of the fluid 

rotational angular velocity 

dimensionless radius (See Chandrasekhar 1939) 

dimensionless density derivative (See Chandrasekhar 1939) 

Inertia tensor 

Kronecker Symbol 



radius of the unperturbed spherical polytrope 

= *„(!-§) 
= GM'/r> 
/j R /(7rGp) 
mean density 

dimensionless angular velocity 
oblateness 

quadrupolar gravitational moment 

Love number (half the apsidal motion constant) 

stellar luminosity out of transit 

stellar luminosity variation in transit 

relative transit depth 

transit radius 

non spherical contribution to the transit depth 



1 

5 

6 

11 

12 

16 

17 

3 



27 



5 

5, 7 
6,8 
6,9 
10 

14, 30 
13 

16, 19 

12 

15 

7 

7 

18 

19 

22 
24 

25, 26 
25,26 
25,26 

31 

36 
37 
42 
48 
48 
48 
48 
50 
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Table B.2. Model parameters for non-irradiated planets of various masses. For each planetary mass (M p ) and age, this table gives 
the spherical equilibrium radius (/?o, p , R] = 7.1492 x 10 7 m), the effective temperature ( r e ff), the central pressure ( P c ), temperature 
( T c ) and density ( p c ) along with the polytropic index (n), the dimensionless moment of inertia (k„, see Eq. (30)) and Love number 
fa). It2 was derived both with our analytical formula (Eq. (42); analytic) and with a numerical method (See Appendix A; numeric). 
An electronic version of this table is available at http://perso.ens-lyon.fr/jeremy.leconte/JLSite/JLsite/Exoplanets_Sitnulations.html 



M p (Mj) Age (Gyr) R , p (R s ) T^jK) P c (MBar) TAK) p c (kg/m 3 ) n K n k 2 

analytic numeric 



0.35 


0.05 


1.1579 


235.6 


2.87 


15 324. 


1016. 


0.913 


0.681 


0.568 


0.428 


0.35 


0.10 


1.0923 


194.6 


3.26 


13 566. 


1 109. 


0.868 


0.694 


0.598 


0.478 


0.35 


0.50 


0.9934 


118.4 


4.05 


9 698. 


1292. 


0.801 


0.715 


0.645 


0.578 


0.35 


1.00 


0.9801 


93.4 


4.18 


9142. 


1326. 


0.792 


0.718 


0.651 


0.594 


0.35 


5.00 


0.9501 


58.9 


4.50 


7 758. 


1402. 


0.771 


0.725 


0.667 


0.631 


0.5 


0.05 


1.1730 


280.9 


5.40 


19045. 


1394. 


0.946 


0.670 


0.547 


0.445 


0.5 


0.10 


1.1112 


234.8 


6.06 


16 765. 


1507. 


0.897 


0.686 


0.579 


0.495 


0.5 


0.50 


1.0165 


140.2 


7.47 


12 120. 


1742. 


0.838 


0.704 


0.618 


0.578 


0.5 


1.00 


0.9961 


108.4 


7.91 


10 950. 


1811. 


0.834 


0.705 


0.621 


0.595 


0.5 


5.00 


0.9682 


62.4 


8.53 


9160. 


1918. 


0.824 


0.708 


0.628 


0.624 


0.6 


0.05 


1.1809 


302.0 


7.67 


21971. 


1670. 


0.987 


0.658 


0.522 


0.443 


0.6 


0.10 


1.1233 


254.0 


8.51 


19 377. 


1795. 


0.936 


0.674 


0.553 


0.490 


0.6 


0.50 


1.0292 


153.9 


10.44 


13816. 


2067. 


0.869 


0.694 


0.597 


0.571 


0.6 


1.00 


1.0063 


120.0 


11.06 


12 332. 


2146. 


0.861 


0.697 


0.602 


0.591 


0.6 


5.00 


0.9775 


66.7 


12.01 


10196. 


2278. 


0.856 


0.698 


0.606 


0.616 


0.7 


0.05 


1.1842 


322.0 


10.44 


25 065. 


1945. 


1.021 


0.648 


0.501 


0.440 


0.7 


0.10 


1.1313 


271.9 


11.50 


22040. 


2090. 


0.974 


0.662 


0.529 


0.484 


0.7 


0.50 


1.0390 


166.9 


14.00 


15781. 


2 398. 


0.901 


0.684 


0.576 


0.562 


0.7 


1.00 


1.0155 


129.9 


14.79 


13 766. 


2492. 


0.887 


0.688 


0.585 


0.583 


0.7 


5.00 


0.9849 


71.0 


16.13 


11251. 


2 644. 


0.882 


0.690 


0.589 


0.607 


0.8 


0.05 


1.1857 


339.9 


13.64 


28 043. 


2218. 


1.044 


0.641 


0.488 


0.439 


0.8 


0.10 


1.1360 


285.7 


14.98 


24 692. 


2 374. 


1.003 


0.654 


0.512 


0.478 


0.8 


0.50 


1.0466 


178.8 


18.13 


17 626. 


2728. 


0.929 


0.676 


0.558 


0.553 


0.8 


1.00 


1.0224 


139.5 


19.22 


15 272. 


2 846. 


0.912 


0.681 


0.568 


0.574 


0.8 


5.00 


0.9903 


76.1 


20.96 


12 243. 


3 022. 


0.903 


0.684 


0.575 


0.599 


0.9 


0.05 


1.1864 


355.5 


17.22 
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1.057 


0.637 


0.480 


0.438 


0.9 


0.10 


1.1393 


297.1 


18.84 


27 295. 


2641. 


1.021 


0.648 


0.501 


0.475 


0.9 


0.50 


1.0525 


189.6 


22.83 


19 399. 


3 048. 


0.952 


0.669 


0.543 


0.546 


0.9 


1.00 


1.0272 


148.9 


24.30 


16774. 


3191. 


0.933 


0.674 


0.555 


0.567 


0.9 


5.00 


0.9946 


82.0 


26.53 


13 201. 


3 402. 


0.921 


0.678 


0.563 


0.592 


1.0 


0.05 


1.1873 


369.1 


21.02 


33 873. 


2 696. 


1.061 


0.636 


0.478 


0.440 


1.0 


0.10 


1.1419 


306.5 


23.04 


29 826. 


2 897. 


1.030 


0.645 


0.496 


0.474 


1.0 


0.50 


1.0576 


199.5 


27.98 


21 165. 


3 357. 


0.970 


0.663 


0.532 


0.540 


1.0 


1.00 


1.0316 


157.6 


29.88 


18218. 


3 528. 


0.951 


0.669 


0.544 


0.560 
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5.00 


0.9980 


87.0 


32.77 


14132. 


3 779. 


0.936 


0.673 


0.553 


0.585 
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0.05 


1.1902 


397.7 


29.16 


39189. 


3127. 


1.051 


0.639 


0.484 


0.445 


1.2 


0.10 


1.1467 


328.1 


32.22 


34 685. 


3 384. 


1.032 


0.645 


0.495 


0.475 


1.2 


0.50 


1.0666 


216.9 


39.33 


24 692. 


3 949. 


0.992 


0.657 


0.519 


0.530 


1.2 


1.00 


1.0396 


173.3 


42.28 


21047. 


4173. 


0.975 


0.662 


0.529 


0.550 


1.2 


5.00 


1.0037 


96.0 


46.87 


15917. 


4513. 


0.960 


0.666 


0.538 


0.574 


1.5 


0.05 


1.1974 


437.2 


43.27 


46493. 


3761. 


1.019 


0.648 


0.502 


0.456 


1.5 


0.10 


1.1541 


362.0 


48.34 


41 254. 


4096. 


1.016 


0.650 


0.504 


0.479 


1.5 


0.50 


1.0768 


238.7 


59.61 


29797. 


4811. 


0.998 


0.655 


0.515 


0.524 


1.5 


1.00 


1.0496 


193.9 


64.49 


25 342. 


5 108. 


0.989 


0.658 


0.520 


0.541 


1.5 


5.00 


1.0103 


108.2 


72.65 


18 523. 


5 584. 


0.976 


0.661 


0.528 


0.563 


1.8 


0.05 


1.2044 


467.7 


59.79 


53 447. 


4405. 


0.994 


0.656 


0.517 


0.465 


1.8 


0.10 


1.1613 


392.9 


67.23 


47 612. 


4816. 


0.999 


0.655 


0.514 


0.484 


1.8 


0.50 


1.0853 


257.2 


83.49 


34 881. 


5 679. 


0.997 


0.655 


0.515 


0.521 


1.8 


1.00 


1.0575 


211.6 


90.76 


29 639. 


6049. 


0.993 
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0.518 


0.535 


1.8 


5.00 


1.0157 


119.2 


103.42 
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0.984 


0.659 


0.524 


0.556 


2.1 


0.05 
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0.980 


0.660 


0.526 


0.471 
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0.10 
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0.50 
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121.62 
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0.530 
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0.05 
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1.1803 
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3.0 
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0.652 


0.508 
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3.0 


5.00 
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290.80 


32425. 


11396. 
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0.652 


0.509 


0.534 
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Table B.3. End of Table B.2. 



Leconte et al.: Distorted extrasolar planets 
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325 053. 
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0.373 
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0.05 
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113 070. 
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Table B.4. End of Table B.2. 



Leconte et al.: Distorted extrasolar planets 



A/t I A /t \ 

M p (Mj) 


Age (Gyr) 


^0,p Cftj) 


^eff (K) 


r c (M car) 


1 C {K) 


Pc (Kg/m ) 


n 




/< 

analytic 


2 

numeric 


25. 


0.05 


1.5200 


2127. 
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815455. 
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0.570 


0.362 


0.395 
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0.10 
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0.50 
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0.565 
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85 812. 
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0.546 


0.325 
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45. 


0.50 
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0.544 


0.322 


0.356 


45. 
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45. 
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797. 
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0.539 
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0.05 
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1.353 


0.552 


0.333 


0.372 


50. 


0.10 
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2443. 
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121580. 


1.382 


0.544 


0.321 


0.352 


50. 


0.50 
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1657 050. 


304 354. 


1.394 


0.540 


0.316 


0.350 


50. 


1.00 


0.9212 


1416. 


201021. 


1 473 950. 


405 043. 


1.401 


0.538 


0.313 


0.350 


50. 


5.00 


0.8133 


877. 


337 922. 


939 025. 
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1.420 


0.533 


0.305 


0.347 


55. 


0.05 
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2 692. 


27 233. 


1 961 060. 


85 721. 
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0.549 
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55. 


0.10 
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44 642. 
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0.542 
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0.348 


55. 


0.50 
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0.537 
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55. 
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0.9152 
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2453 178. 
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0.533 


0.305 
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1.00 


0.9040 
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0.537 


0.310 
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Leconte et al.: Distorted extrasolar planets 

Table B.5. Same as TableB.2, for the most irradiated planets {F+ = 4.18 X 10 6 W.rrT 2 ). An electronic version of this table is 
available at http://perso.ens4yon.fr/jeremy.leconte/JLSite/JLsite/Exoplanets_Simulations.html 
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Age (Gyr) 
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0.9 


5.00 


1.1048 


96.5 


21.02 


22 830. 


2 866. 


0.981 


0.660 


595 

U. Ji J 


479 

U.'T / y 


1.0 


0.05 


1.3845 


372.0 


14.45 


45 609. 


2 005. 


1.155 


0.608 


0.427 


0.322 


1.0 


0.10 


1.3074 


310.6 


16.76 


42025. 


2248. 


1.124 


0.617 


443 

U.^t^tJ 


357 


1.0 


0.50 


1.1918 


195.7 


21.29 


33 727. 


2719. 


1.060 


0.636 


0.478 


0.424 


1.0 


1.00 


1.1616 


159.4 


22.69 


30526. 


2 863. 


1.035 


0.644 


493 


445 


1.0 


5.00 


1.1038 


100.7 


25.92 


24599. 


3 168. 


0.994 


0.656 


0.517 


0.482 


1.2 


0.05 


1.3802 


403.7 


19.45 


51893. 


2287. 


1.082 


0.630 


466 


347 


1.2 


0.10 


1.3038 


337.8 


22.81 


47 931. 


2582. 


1.070 


0.633 


473 

U.'t 1 J 


379 
u.j / y 


1.2 


0.50 


1.1886 


213.2 


29.76 


38 662. 


3176. 


1.048 


0.640 


485 


437 

U.T - J / 


1.2 


1.00 


1.1592 


173.3 


31.97 


35 060. 


3 362. 


1.034 


0.644 


494 


454 

U.T - J^t 


1.2 


5.00 


1.1032 


108.7 


36.81 


28 050. 


3 754. 


1.007 


0.652 


510 

u — uu 


487 


1.5 


0.05 


1.3816 


441.7 


28.32 


59 797. 


2739. 


1.010 


0.651 


508 

U.JuO 


375 

U. J / J 


1.5 


0.10 


1.3063 


371.5 


33.41 


55 635. 


3 094. 


1.009 


0.652 


508 

U. JUO 


403 


1.5 


0.50 


1.1899 


236.2 


44.60 


45 509. 


3 844. 


1.017 


0.649 


503 

u. juj 


453 


1.5 


1.00 


1.1602 


191.8 


48.41 


41202. 


4100. 


1.016 


0.649 


504 
u. jut 


466 


1.5 


5.00 


1.1047 


119.5 


56.43 


32 981. 


4612. 


1.005 


0.653 


511 

U.J 1 1 


494 

U.t74 


1.8 


0.05 


1.3842 


476.2 


39.02 


67811. 


3 207. 


0.976 


0.662 


578 

U.JiO 


394 
u. jyi 


1.8 


0.10 


1.3099 


401.5 


46.16 


62 943. 


3 636. 


0.976 


0.661 


578 

U.JiO 


0.421 


1.8 


0.50 


1.1931 


257.0 


62.12 


51794. 


4534. 


0.993 


0.656 


518 

u — no 


464 

U.T-U^t 


1.8 


1.00 


1.1625 


208.6 


67.92 


47 144. 


4853. 


1.000 


0.654 


0.514 


0.476 


1.8 


5.00 


1.1068 


129.5 


79.87 


37 666. 


5 492. 


0.999 


0.655 


0.514 


0.499 


2.1 


0.05 


1.3877 


507.3 


51.29 


77 544. 


3 643. 


0.963 


0.665 


0.536 


0.408 


2.1 


0.10 


1.3137 


429.3 


61.02 


71 611. 


4 159. 


0.964 


0.665 


0.536 


0.432 


2.1 


0.50 


1.1963 


276.7 


82.88 


58167. 


5 246. 


0.980 


0.660 


0.526 


0.472 


2.1 


1.00 


1.1650 


224.2 


90.91 


53 093. 


5 622. 


0.990 


0.657 


0.520 


0.482 


2.1 


5.00 


1.1087 


138.8 


107.73 


42 307. 


6 398. 


0.996 


0.656 


0.516 


0.502 


3.0 


0.05 


1.3997 


583.5 


96.63 


106857. 


4903. 


0.960 


0.666 


0.538 


0.435 


3.0 


0.10 


1.3263 


499.2 


116.54 


99 725. 


5 658. 


0.964 


0.665 


0.536 


0.453 


3.0 


0.50 


1.2058 


329.1 


164.27 


80267. 


7 377. 


0.983 


0.659 


0.524 


0.480 


3.0 


1.00 


1.1709 


267.7 


183.21 


72 047. 


8 026. 


0.993 


0.657 


0.518 


0.488 


3.0 


5.00 


1.1108 


163.4 


222.49 


56621. 


9310. 


1.006 


0.653 


0.510 


0.501 
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Table B.6. End of Table B.5. 



Leconte et al.: Distorted extrasolar planets 



Mp(Mj) Age (Gyr) i? , P P c (MBar) TAK) p c (kg/m J ) n~ * fc 2 

analytic numeric 



S 


n n5 


1 .4224 


741.1 


247.80 


166 0m 

1 UU 71U. 


7 87n 
/ o / u. 


n 007 

U.77 / 


n 655 


n 5 5 

U.J 1 J 


n 457 




0.10 


1 1471 


610 5 


in4 77 

JUT - . / / 


157 781. 


??0 
y zz;?. 


1 nn7 

1 .UU / 


n 65? 


n 5no 

U.JU7 


466 


S 


n 5n 


1 ?1 13 




464 41 


1 ?8 666 

1 ZO UUU. 


12786. 


1.044 


n 641 

U.U^t 1 


n 488 


n 476 
u.^t / u 


S 


i on 


1 1 7ns 

1.1/ uo 


344.7 


517 ?4 


1 15 271. 


14?01 


1 n50 

1 .UJ7 


0.637 


n 470 
u.t - / y 


n 478 

U.^T / O 


S 

J.VJ 




1.1014 


7118 


684.92 


on 171 

y\j i / 1 . 


17 174. 


1.077 


n 611 

U.U J 1 


n 460 


n 481 

U.'tO J 


8 

O.V/ 


n n5 


1 .4420 


01? 1 

7Ji.l 


616 18 


?50 574 


1 ? 041 

1 Z 7tj. 


1 ns6 

1 .uou 


0.628 


n 464 


n 440 


O.V/ 


n in 


1 1405 


701 n 


811 41 

O J 1 .T - J 


?48 1?6 


1 5 866 
1J ouu. 


1 m5 

1 . 1UJ 


n 6?i 


n 451 


n 451 


O.V/ 


n 5n 

U. JU 


1 ?ns 1 

1 .iUJ 1 




1 327.72 


212 307. 


?? 576 

ZZ J / U. 


1 140 

1 . It7 


n 6in 

U.U 1U 


n 4in 


n 451 


8 

O.V/ 


1.00 


1.1571 


445 5 


1 57n 48 

1 J / U.T-O 


192724. 


?5 501 

zj jyj. 


1.164 


n 6n5 


0.422 


n 451 

U.H-J 1 


8 

O.V/ 


s nn 


1 .U / T-O 


278.7 


2 1 1 1.44 


147 177. 


11 871 

J 1 O / J. 


1 1 70 

1.1/7 


n 6ni 

U.UU 1 


0.414 


n 45? 




n n5 


1 44.81 


1 OS 1 s 

1 UJ 1 .o 


1 016.77 


1?5 1 64 


1 6 611 

1 U U J J. 


1.143 


n 6i? 

U.U li 


n 411 


n 418 


1 


n in 


1.3480 


801 7 


1 359.72 


111118 

J 1 J J 1 o. 


20717. 


1 1 61 

1 . 1 U J 


n 6n6 

U.UUU 


0.422 


n 44n 


10 


n 5n 

U. JU 


1 1 01 1 

1 . 1 7J 1 


588 7 


2 240.54 


?7n4on 

Z / U t7U. 


inn70 

ju u / y. 


1 ?nn 


n 505 


n 4ni 


n 418 


10 


i nn 

1 .UU 


1 i4sn 


5nn 6 


2 647.31 


247482. 


14ni1 


1.210 


n 50? 


n ioo 


n 417 


1 


s nn 


1 n58i 

1 .UJOJ 


321.2 


1616 15 

J U1U. 1 J 


1 88 1 00 

1 OO 1 77. 


42 848. 


1.217 


n 5on 

U.J7U 


n 105 


n 417 


1 3 


n n5 


1 4541 


1 227.3 


1 771.70 


4?7 6?6 
H'Z / uzu. 


22 312. 


1.207 


n 501 


n 4nn 


0.422 


1 3 


0.10 


1 i4?n 


1 036.9 


2453.77 


414 870. 


28 462. 


1.228 


n 587 

U.JO / 


n ion 


0.422 


1 3 


n 5n 


1.1746 


684 1 

UO't. 1 


4 190.16 


16n 16? 


4? Ill 
^tz jjj. 


1.249 


n 58i 

U.JO 1 


n i8n 

U. JOU 


0.421 


1 3 


i nn 


1 1 ?fin 

1 . 1 Z.UU 


57? 1 

J / z.. 1 


4 05? 80 


in no5 

Jjl U7J . 


47 onn 

t / yuu. 


1 ?5 1 

1 .Z J 1 


n 58n 


n 170 

u. J / y 


n 4? 1 

U.H-Z 1 


13.0 


5.00 


1.0376 


378.2 


6 812.26 


253 286. 


60 644. 


1.245 


0.582 


0.381 


0.422 


17.0 


0.05 


1.9236 


1814.9 


1 092.21 


530529. 


13 540. 


1.206 


0.593 


0.401 


0.392 


17.0 


0.10 


1.5737 


1 547.6 


2 343.85 


565 858. 


24 210. 


1.250 


0.581 


0.379 


0.402 


17.0 


0.50 


1.1635 


839.0 


7 729.08 


492 879. 


58 757. 


1.285 


0.571 


0.363 


0.405 


17.0 


1.00 


1.1059 


676.3 


9432.93 


451 135. 


68 127. 


1.280 


0.572 


0.365 


0.406 


17.0 


5.00 


1.0153 


451.7 


13 139.79 


350018. 


87 181. 


1.267 


0.576 


0.371 


0.409 


20.0 


0.05 


1.7413 


1 796.7 


2 127.33 


643 444. 


20 955. 


1.257 


0.579 


0.376 


0.392 


20.0 


0.10 


1.4005 


1 444.2 


5 080.99 


666543. 


40 165. 


1.300 


0.566 


0.356 


0.395 


20.0 


0.50 


1.1447 


901.3 


11241.96 


577 707. 


72 520. 


1.299 


0.567 


0.357 


0.398 


20.0 


1.00 


1.0891 


731.1 


13 657.81 


530282. 


83 878. 


1.292 


0.569 


0.360 


0.399 


20.0 


5.00 


1.0015 


492.3 


18 915.42 


418 729. 


106 994. 


1.277 


0.573 


0.366 


0.402 



